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single  crystal  of  CdIn2Te^  was  grown  near  the  end  of  the  quarter.  In 
addition,  synthesis  of  ZnIn2S4  was  started  using  the  constituent  Ele¬ 
ments  but  failed  to  go  to  completion  because  a  sulfide  layer  formed  on 
top  of  the  molten  metals  and  prevented  further  reaction  with  gaseous 
sulfur. 
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Quantitative  calculations  of  EO  coefficients  were  made  for  ternary  com¬ 
pounds  and  showed  close  agreement  with  measured  results. 
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REPORT  SUMMARY 


The  objectives  of  this  program  are  to  find  and  develop  new  IR 
transmitting  materials  and  to  provide  new  data  on  the  electrooptic  (EO) 
properties  of  those  most  likely  to  have  an  EO  coefficient  an  order  of 
magnitude  higher  than  materials  currently  in  development  for  tunable 
filters.  The  mai-  technical  problems  anticipated  include  the  synthesis 
and  single-crystal  growth  of  these  materials:  many  are  poorly  charac¬ 
terized  and  others  have  high  melting  points  or  melt  incongruently .  Our 
approach  will  overcome  these  obstacles  by  first  synthesizing  20  poly¬ 
crystalline  samples;  subsequently,  dielectric  constants  at  low  and 
ambient  temperatures  will  be  determined,  and  the  two  best  materials  of 
the  survey  will  be  grown  as  single  crystals  (second  year  of  the  program). 

During  the  last  quarter,  emphasis  was  placed  on  the  growth  of 
single-crystal  Cdl^Te^.  Using  the  published  phase  equilibrium  diagram, 
the  liquid-solid  relationships  were  established  for  crystal  growth.  A 
small  single  crystal  of  CdIn2Te^  was  grown  near  the  end  of  the  quarter. 

In  addition,  synthesis  of  Znl^S^  was  started  using  the  constituent 
elements  but  failed  to  go  to  completion  because  a  sulfide  layer  formed 
on  top  of  the  molten  metals  and  prevented  further  reaction  with  gaseous 
sulfur. 

Quantitative  calculations  of  EO  coefficients,  based  on  the  bond 
charge  dielectric  theory  of  Phillips  and  Van  Vechten,  were  made  for 
ternary  compounds  and  showed  good  agreement  with  measured  results. 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 

A.  PROGRAM  OBJECTIVES 

The  objectives  of  this  program  are  to  find  and  develop  new  IR 
transmitting  materials  and  to  provide  new  data  on  the  electrooptic  (EO) 
properties  of  those  most  likely  to  have  EO  coefficients  an  order  of 
magnitude  higher  than  materials  currently  in  development  for  tunable 
filters.  The  main  technical  problems  anticipated  include  the  synthesis 
and  single-crystal  growth  of  these  materials:  many  are  poorly  charac¬ 
terized  and  others  have  high-melting  points  or  melt  incongruently .  Our 
approach  will  overcome  these  obstacles.  First,  we  will  synthesize  20 
polycrystalline  samples.  Then  the  dielectric  constant  of  each,  at  both 
low  and  ambient  temperatures,  will  be  determined,  and  the  two  best  mate¬ 
rials  of  the  survey  will  be  grown  as  single  crystals  (second  year  of 
the  program) . 

B.  SUMMARY 

During  the  last  quarter,  emphasis  was  placed  on  the  growth  of  single¬ 
crystal  Cdln^Te^.  Using  the  published  phase  equilibrium  diagram,  the 
liquid-solid  relationships  were  established  for  crystal  growth.  A  small 
single  crystal  of  CdIn2Te^  was  grown  near  the  end  of  the  quarter.  In 
addition,  synthesis  of  ZnIn?S^  was  started  using  the  consituent  elements 
but  failed  to  go  to  completion  because  a  sulfide  layer  formed  on  top  of 
the  molten  metals  and  prevented  further  reaction  with  gaseous  sulfur. 

Quantitative  calculations  of  EO  coefficients,  based  on  the  bond 
charge  dielectric  theory  of  Phillips  and  Van  Vechten,  were  made  for 
ternary  compounds  and  showed  good  agreement  with  measured  results. 
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SECTION  2 


MATERIALS  PREPARATION  AND  CRYSTAL  GROWTH 

A.  Cdli^Te^ 

Primary  emphasis  during  this  quarter  was  placed  on  the  growth  of 

single-crystal  CdIn2Te^,  a  compound  that  we  previously  had  reported^  as 

having  a  high  value  for  the  low-frequency  dielectric  constant.  Using 

the  solid-liquid  relationships  shown  (Figure  1)  in  the  pseudo-binary 

2 

phase,  diagram  CdTe-I^Te^  as  published  in  the  literature,  we  attempted 
to  grow  single-crystal  material.  Melt  compositions  were  selected  lying 
in  the  composition  range  covered  by  the  line  a-b.  In  this  region,  solid 
Cdl^Te^  (3  =  solid  Cdl^Te^)  is  in  equilibrium  with  liquid  in  the  tem¬ 
perature  range  785°C  to  702°C.  By  selecting  a  starting  compositon  close 
to  "a,"  crystals  of  pure  Cdl^Te^  can  be  grown  by  lowering  the  tempera¬ 
ture.  Below  702°C,  however,  additional  phases  begin  to  solidify  and 
cause  a  multiphase  polycrystalline  region  to  begin  to  grow,  thus  pre¬ 
venting  further  single-crystal  growth. 

Several  runs  were  made  utilizing  the  conditions  discussed  above 
to  grow  single-crystal  CdIn2Te^.  Fused  quartz  crucibles  approximately 
16  mm  in  diameter  with  conical  bottoms  were  used  to  nucleate  single¬ 
crystal  growth.  The  sealed,  evacuated  ampoules  were  processed  in  a 
vertical  two-zone  Bridgman-like  furnace.  The  ampoule  was  rotated  and 
vibrational  stirring  was  applied  during  the  growth  cycle.  As  antici¬ 
pated  from  the  above  discussion,  several  ingots  were  produced  in  which 
the  lower  (conical-shaped)  sections  were  single  phase  and,  in  one  case, 
a  single  crystal  of  CdIn2Te^;  the  upper  parts  were  multiphase  and 
polycrystalline . 

X-ray  microprobe  analysis  of  the  single-phase  material  showed 
insignificant  variations  in  composition  from  nose  to  tail  with  an 
average  composition: 

58%  Te 

29%  In 

13%  Cd. 
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Figure  1.  Phase  equilibrium  diagram  for  CdTe-In^Te^. 


There  are  indications  from  the  analysis  that  the  composition  is 
slightly  rich  in  In^Te^  and  that  there  is  some  solid  solubility  between 
Cdln^Te^  and  In^e^.  However,  our  results  indicate  that  the  solid  solu¬ 
bility  line  between  the  regions  is  more  vertical  than  that  reported  in 
Ref.  1  (dashed  line  in  Figure  1)  and  probably  closer  to  the  vertical 
(dotted)  line  in  Figure  1. 

One  of  the  most  successful  ingots  from  these  runs  (Figure  2) 
yielded  a  single-crystal  section  (cone-shaped)  as  well  as  the  final 
multiphase  section;  the  boundary  between  the  sections  is  noted  in  Fig¬ 
ure  2.  The  tip  of  the  cone  appears  as-grown  and  indicates  a  poor  start 
with  nucleation  getting  underway  on  an  apparent  gas  void,  resulting  in 
a  concave  interface  observed  during  crystallization.  That  the  crystal 
grew  single  under  these  conditions  (see  Figure  3,  SEM  photo  of  as-grown 
tip),  as  confirmed  by  the  growth  patterns  shown  in  Figure  3  and  by  Laue 
patterns,  is  very  encouraging  for  future  single-crystal  growth  runs. 

Plans  for  the  next  quarter  include  further  attempts  at  single¬ 
crystal  growth  of  CdIn2Te2  seeking  optimum  temperature- lowering  rates 
to  prevent  the  crystal  from  cracking  as  well  as  solution  from  being 
entrapped  to  obtain  the  largest  high-quality  single  crystal  for  evaluation 
of  its  RO  properties. 

B.  ZnIn„S. 

2  4 

During  the  last  quarter,  considerable  effort  was  expended  to  syn¬ 
thesize  Znl^S^  from  its  constituent  elements:  Zn,  In,  and  S.  A  long¬ 
term  run  in  which  the  molten  metals  were  exposed  to  gaseous  sulfur  failed 
to  go  to  completion  after  several  weeks  and  was  aborted.  Subsequent 
observation  indicated  that  a  layer  of  the  sulfide  had  formed  on  top  of 
the  molten  metals  and  prevented  further  reaction.  A  second  run  has 
gotten  underway  with  added  iodine  to  initiate  vapor  transport  of  the 
solidified  species  and  allow  further  reaction  to  take  place. 


Figure  2.  Ingot  containing  single-crystal  CdIn„Te,. 

(Region  to  right  of  dashed  line  is  single.) 
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MATERIALS  EVALUATION 
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A.  PROPERTIES  OF  Cdln-Te, 

l  4 

3-7 

Cdli^Te^  has  been  reported 

•  Crystal  class 

•  Band  gap 

•  Resistivity  (r.t.) 

•  Mobility 


to  have  the  following  properties: 
Tetragonal,  4 
0.9  eV 

Maximum  10"*  ft-cm  (n  type) 

Minimum  60  £2-cm 

4000  cm^V  '''sec  3 


•  Transmission  range  1  to  37  Mm 

•  Density  5.9  g/cm3. 

Our  measurements  indicate  that  the  sample  we  grew  has  a  resistivity 
of  107  ft-cm  at  room  temperature.  The  dielect  ic  constant  at  10  kHz  was 
measured  to  be  456  at  room  temperature.  Low- temperature  measurements 
showed  a  significant  decrease  in  dielectric  constant  to  73  at  83°K.  It 
is  possible  that  Cdln^Te^  is  a  ferroelectric  slightly  below  room  tem¬ 
perature.  Our  evaluation  to  date  has  not  proven  that  it  is,  but  there 
is  no  definite  evidence  that  it  is  not.  This  evaluation  is  continuing, 
and  our  results  will  be  discussed  in  a  future  report. 


B.  QUANTITATIVE  CALCULATION  OF  EO  COEFFICIENTS 
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The  following  represents  the  work  of  Professor  Amnon  Yariv  and  his 
graduate  student,  C.  Shih,  of  the  California  Institute  of  Technology.  They 
have  applied  the  bond  charge  dielectric  theory  of  Phillips  and  Van  Vechten 
to  the  calculation  of  the  E0  tensor  coefficients.  Comparisons  with 
experimental  values  for  binary  compounds  having  zincblende  and  wurtzite 
structures  (as  previously  reported,  Quarterly  Report  4)  is  very  good. 
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Yariv  and  Shih  recently  extended  the  calculations  to  ternary  compounds 
e.g.,  LiNbO^,  CdIn2Te^,  where  they  also  find  very  good  agreement.  These 
results  will  appear  in  Table  I.  A  discussion  of  the  calculation  appears 
in  the  appendix. 

We  applied  the  theoretical  model  for  calculating  the  electrooptic 
coefficient  to  the  crystal  CdIn2Te^.  We  used  the  following  data: 

O 

•  Covalent  radii  (Cd.In.Te)  —  1.405A 

•  Structure  —  4 

Using  published  data  of  bond  susceptibilities,  we  calculated  the 
optical  dielectric  constant  =  7.25.  Using  the  above  data  in  the 
equations 


(r .  )  . 

ljk  ionic 


£o(Cdck  "  c”k)  y  _6 

V  Ne*  c'  c' 
c  i  j 


+  l/2(<6jk  + 

=  -  1 .As) f ^  -  0.02  , 


(1) 


we  obtain  the  result 


r. ,  =  220  x  10 
41 


-12 


m/V  ±30% 


This  predicted  value  is  a  huge  number  (roughly  100X  that  of  GaAs  and 
that  of  LiNbO^).  The  main  differ< 
as  an  example,  is  due  to  the  factor: 


10X  that  of  LiNbO^).  The  main  difference  between  this  crystal  and  LiNbOy 


dc 


EX  , 

l  .) 
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Table  I.  Results  for  LiNbO^  and  LiTaO^ 


LiNbO 

3 

LiTa03 

'  dc3 

28 

43 

rdcl,2 

43 

41 

e  C/E 

1 

.8 

2.0 

Nb-0  (short) 

Nb-0  (long) 

Ta-0  (short) 

Ta-0  (long) 

2r 

o 

1.889  A 

2.112 

1.891 

2.071 

fi 

0.821 

0.830 

0.847 

0.853 

-f 

0.292 

0.241 

0.282 

0.238 

r 

33 

r51 

r33 

r51 

r . 

ionic 

+19.9 

+19.7 

+27.8 

+16.5 

r  - 
elec 

+  6.0 

+  0.8 

+  3.70 

+  0.24 

theo 

r 

sum 

+25.9 

+20.5 

+31.5 

+16.7 

r 

exptl 

+28 

+23 

+30 

+15 

which  is  equal  to  8.55  in  CdIn„Te.  and  to  1.02  in  LiNbO_.  It  will  be 

2  4  3 

very  exciting  to  see  if  this  prediction  is  actually  true. 

We  also  predict  a  large  EO  coefficient  for  Znln,-,S^  since  it  has 
3m  symmetry  so  that  the  form 


I- 

u 


is  non-vanishing  for  r..,.  The  more  numerical  orediction  awaits  the 
result  of  measurement  of  e<jc* 
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APPENDIX 

QUANTITATIVE  CALCULATION  OF  ELECTROOPTIC  COEFFICIENTS 


The  second-order  nonlinear  optical  response  of  asymmetric  crystals 
is  usually  represented  by  the  relation  pf  +  ^  ^  ^  between  the 

amplitude  of  the  induced  polarization  at  to  +  ft  and  the  inducing  field 
amplitudes  at  to  and  ft.  The  case  when  both  to  and  i  are  optical  frequen¬ 
cies  (i.e.,  they  are  frequencies  above  the  lattice  response  but  below 
the  optical  absorption)  has  been  considered  by  Levine.  He  used  the 
localized  bond  charge  model  (see  Figure  A-l(a))  of  Phillips  and 
Van  Vechten2  (PV) ,  which  attributes  the  dielectric  response  of  covalent 
crystals  to  the  localized  bond  charge  resulting  in  a  linear 
susceptibility: 


X  = 


(1) 


where  ft  is  the  plasr  frequency  due  to  valence  electrons,  and  E  (the 
p  9  9  2  B 


effective  energy  gap)  is  given  by  E^  =  E^  +  C  (where  E^  is  the  homopolar 
component,  and  C  is  the  heteropolar  (ionic)  component  of  the  gap  energy). 


(oj)  _(to)  ,  . 

Levine  starts  with  the  linear  dielectric  response  P.  -  X,- »  taking 

1  (ft}  j 


v  .  to  be  an  instantaneous  function  of  the  second  field  E, 
ij  K 


This  field 


J 

causes  a  change  Ar  (=  -Arp)  in  the  bond  charge  position,  as  shown  in 
Figure  A-l(a),  that  oscillates  at  ft.  The  explicit  dependence  of  E  and 

9  6 

C  on  r  given  by  PV  is  then  used  to  obtain  X^j^’  where 


-  x'j0)  +  4xt.(t)  -  X<f  +  2X1Jk4‘ ’  ac 


(2) 
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Figure  A-l.  The  response  of  ion  and  bond-change  to  the  applied  elec¬ 
tric  field.  In  (a)  and  (b)  the  frequency  of  the  field 
is  higher  than  the  lattice  response.  Only  the  displace¬ 
ment  of  the  bond  charge  takes  place.  In  (c)  and  (d)  the 
frequency  of  the  field  is  lower  than  the  lattice  response. 
It  induces  the  displacement  of  the  bond  charge  (Ar  ) , 
bond  length  elongation  (Ad^),  and  bond  rotation  (A0T. 
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If  the  frequency  ft  is  below  the  lattice  response  region  (we  will 
refer  to  it  in  this  case  as  "low"),  then,  in  addition  to  the  purely 
electronic  nonlinear  response  described  above,  there  is  now  a  contri¬ 
bution  to  Xi-jk  because  the  crystalline  ions  are  capable  of  following 
the  field  cos  ftt.  This  is  illustrated  in  Figure  A-l(c  and  d).  In 

addition  to  the  displacement  Ar„(  of  the  covalent  bond  charge,  there  is 
now  an  elongation  AdQ  of  the  atomic  separation  and  a  rotation  0  of  the 
bond  direction;  these  are  caused  by  the  ionic  displacement  Axk-  Ax^ 
is  obtained  from  "low"  frequency  dielectric  constant  measurements  and 

is  used  to  obtain  Ad  and  A9.  We  use,  in  the  spirit  of  Levine,1  the 

o 

change  Ad  to  calculate  the  corresponding  change  Ay  cos  ftt  in  bond  sus- 
o 

ceptibility.  This  will  give  rise  to  a  polarization 

(artft)  _  p(w)  (ft) 

Pi  ‘  Xijkbj  Ek 


A  second  contribution  to  X- -u  is  due  t0  the  rocking  at  ft  of  the  bond 

*•*  3 

angle  (0  =  0q  +  A0  cos  ftt),  which  yields  a  dipole  component  along  i  at 
(w4ft)  even  when  AdQ  =  0. 

In  what  follows,  we  will  obtain  expressions  for  the  ionic  contri¬ 
bution  to  X..,  »  which  is  due  to  Ad  and  A0.  Adding  this  result  alge- 
ijk  o 

braically  to  xf-^Cas  measured  by  second  harmonic  generation  experiments 
i  i  k 

or  calculated  by  Levine  )  yields  the  total  nonlinear  tensor 


ionic  elect 
Xijk  “  Xijk  Xijk 


The  constants  thus  determined  are  those  which  characterize  the 
ijk 

line..r  electrooptic  (Pockels)  effect.  The  relationship  between  the 
conventionally  defined  electrooptic  tensor  r  and  XAjk  is  (in  fflcs 
units) 


:wwun  w  i/v.  %vvMHHBMHMHHMHHHnnMnMMnNsnmnnnnnnni 


e  e 


4jk 


2e 


ijk 


(3) 


We  use  this  procedure  below  to  calculate  the  tensor  r^j  in  a  number  of 


zincblende  and  wurtzite  crystals.  The  results  are  in  good  agreement 
with  experiments. 

The  linear  susceptibility  of  a  diatomic  crystal  is  given  as 

The  ionicity  and  covalency  of  the  bond  are 

respectively.  Values  of  E^,  C,  f^, 

The  expres- 


X  =  (hfl  )2/eI  (Ref.  2) 


P, 


defined  as  f.  =  C^/E^  and  f  =  E^/E^ 
i  g  c  h  g 


and  f.  for  a  large  number  of  crystals  are  given  in  Ref.  3. 


sions  used  in  the  evaluation  of  E,  and  C  are 

n 


1,2 


"h 


s  =  2.48 


(4) 


-k 

C  «  e  s 


Va  r3  / 


(5) 


a 


rS  ro 


-ksro 


where  d  =  r  +r„  is  the  bond  length,  r  „  are  the  atomic  radii,  e 

o  a  g  °  a, 3 

is  the  Thomas-Fermi  screening  factor,  and  the  proportional  factors  in  E^ 
and  C  are  independent  of  bond  length  and  atomic  radius. 


To  consider  the  crystals  with  highly  unequal  atomic  radii,  a 
generalized  form  for  E^  was  proposed:'*' 


-2 


2s 


[<V 


.2s  .  .2s 

+  (r6-rc)  , 


2«2  -r  ) 
o  c 


2s 


(6) 


where  r  is  the  average  core  radius,  which  is  included  here  since  the 


contribution  to  the  susceptibility  in  the  core  region  is  very  small. 
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The  linear  macroscopic  susceptibility  tensor  x^j  is  related  to  the 


bond  polarizability  by 


Xij  vXanian/n 


where  V  is  the  volume  of  the  unit  cell,  a  ,  is  the  direction  cosine  of 
the  bond,  n  refers  to  the  individual  bond,  and  the  summation  is  over  all 
the  bonds  in  a  unit  ceil.  Although  PV  described  the  macroscopic  sus¬ 
ceptibility  X  in  terms  of  the  average  energy  gap,  we  assume  that  li^  and 
C  should  be  related  to  the  bond  polarizability  directly,  i.e., 

(hfi  )2 

6  «  ■  9P  ■  . 
n  E2 
g 

The  change  of  bond  polarizability  due  to  the  applied  electric  field  is 


M^)  -  o 

-f  - +  fcEh  ^  > 


2f  *2 
ri  C 


When  the  bond  length  varies,  it  is  reasonable  to  assume  that  the  ratio 
of  r  to  rQ  remains  constant.  With  this  assumption,  the  two  independent 

OL  p 

parameters,  r  D,  can  be  transformed  into  two  quantities  that  relate 

a,p 

directly  to  the  macroscopic  properties  of  crystals: 


Ar  =  Ad  +6 
ad  o 
o 


Are  =  T  Ad0  *  6  ’ 

o 


where  6  is  the  displacement  of  the  bond  charge  in  the  case  of  no  bond 
elongation  AdQ  =  0.  From  Eqs.  5,  6,  8,  and  9,  the  change  of  bond  polar¬ 
izability  is  obtained  as  (here  we  drop  the  bond  index  n) : 


where  p  =  (r  -rQ)/(r  +ra).  In  the  first  term,  k  r  /2  is  obtained  for 

a  b  d  b  s  o  . 

the  screening  wave  number,  is  proportional  to  d  7  (Ref.  4),  and 

the  number  (-3/2)  is  due  to  the  fact  that  fi  is  proportional  to  d-3^2. 

P  o 

The  second  term  on  the  right  side  of  Eq.  10  is  exactly  the  same  as  the 
expression  obtained  by  Levine  in  his  calculation  of  the  nonlinear  opti¬ 
cal  suscept 'bilitv.5  The  first  term,  which  is  proportional  to  d  ,  is 

o 

thus  the  ionic  contribution  of  a  single  bond  due  to  bond  stretching. 

The  rotational  contribution  can  be  obtained  by  considering  the 
changes  in  bond  direction  cosines.  These  are  related  to  the  ionic  dis¬ 
placement  Ax^  by 


Aa 


ni 


■  (5ik-a„iV> 


From  Eq.  7  we  have 


1  V' 

A X. .  =  - n  (a  a  A8  +  Aa  .a  .8  +  a  ,Aa  .8  ) 

ij  V  n  ni  nj  n  ni  nj  n  nx  nj  n 


(ID 
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The  complete  ionic  contribution  to  the  nonlinear  susceptibility  is  thus 


AX 


ion 

ij 


L  i* 

'  |n  Vro 


fa  .a  4a  , 

ni  nj  nk 


+  2  (ank6jk 


“njV 


Xxk 


(12) 


where 


2.5 


k  r  \ 

-  1A8J  f.  -  0.02  . 


OL3) 


Ax^  is  related  to  the  dielectric  constant  of  the  crystal  as 


ft 

Ne  Ax. 
c  k 


=  e 

o 


(14) 


where  N  is  the  number  of  pairs  of  atoms  per  unit  cell,  e  is  the  Callen 

6  ^ 
effective  ionic  charge,  £■'  is  the  relative  dielectric  constant,  c'  is 

dc 

the  relative  optical  permittivity,  and  is  the  low-frequency  electric 
field  component  along  the  k  direction. 

Using  Eq.  3  we  obtain  the  final  working  expression  for  the  electro¬ 
optic  tensor: 


ion 

'ijk 


e  (e ,  ,  -e  , 

o  dck  °°k 

VNe*£'e' 
c  i  j 


)  (  en  r 

-  '2 —  fa  .a  .a  ,  + 

|n  r0  L  ni  nj  nk 


4  (a  .6  +  a  6  ) 

2  m  jk  nj  ik 


(15) 
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For  wurtzite  crystals,  we  neglect  the  small  distortion  from  the 
perfect  tetragonal  structure.  &n  can  be  expressed  in  terms  of  the  mea¬ 
sured  macroscopic  susceptibility  x  as  in  Eq.  7,  and  the  electrooptic 
coefficients  of  zincblende  and  wurtzite  crystals  are  obtained  as 
follows: 


2  , 
a  wf 

zincblende  r  ,  =  0.3689  — . — 

14  * . 

e  /e 
c 


(16) 


wurtzite  r 


33n  ■  "2  ro"  ■  °-',26°  -T 


wf 


e  /e 
c 


(17) 


where  a  is  the  lattice  constant;  a  cc  =  a2 c  ;  w  =  (e-1)  (e  ,  -e)/e2  ; 

0  eff  oo  dc  ’ 

the  r  are  in  units  of  10  “  m/V;  and  a  ,  a  cc  are  in  units  A.  Values 

J-J  o  *efr 

of  the  parameters  aQ,  a£ff,  w,  f±,  f,  and  e^  are  listed  in  Table  A-I. 

The  dependence  of  the  electrooptic  coefficients  r*°£  on  the  bond 

ij  k 

geometry  is  perhaps  the  most  illuminating  feature  to  emerge  from  this 
work.  This  dependence  is  contained  in  curly  brackets  in  Eq.  15.  For 
diatomic  single-bond  crystals,  B  is  a  constant,  and  the  geometrical 
factor  becomes 


Gijk  =  bolds 


f  a  a  a  +  \  (a  ,  6  +  a  ) 

m  nj  nk  2  nk  jk  nj  ik 


The  factor  f  is  typically  |f|  £  0.3.  Table  A-II  contains  a  listing  of 
these  factors  for  some  key  directions  (ijk)  in  crystals  of  the  zinc¬ 
blende,  wurtzite,  and  LiNbO^  classes. 

It  follows  immediately  that  when  E  a  J  0  the  second  term  in  G 

n  ni  ijk 

is  an  order  of  magnitude  larger  than  the  first  one.  In  such  crystals, 
the  ionic  contribution  to  r  ,  is  about  an  order  of  magnitude  larger 
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Table  A-I.  Parameters  and  Results  of  EO  Calculations 


AB 

Zincblende 

Wurzite 

GaAs 

GaP 

ZnSe 

ZnS 

ZnTe 

CuCl 

ZnS 

CdS 

CdSe 

a 

a 

5. 

65 

5.45 

5.67 

5.41 

6. 

09 

5 

41 

5 

39 

5.85 

6 

08 

Cdc 

13. 

2b 

12. 0C 

9.ld 

8.3d 

10. 

ld 

7 

5e 

8 

7f 

9.4d 

10 

2d 

w 

0. 

192 

0.284 

C.  450 

0.528 

0. 

331 

0 

656 

0 

567 

0.652 

0 

562 

*1* 

0. 

310 

0.370 

0.630 

0.623 

0. 

546 

0 

749 

0 

623 

0.683 

0 

699 

f 

-0. 

091 

-0.113 

-0.163 

-0.179 

-0. 

119 

-0. 

212 

-0 

181 

-0.162 

-0 

147 

e  /e 
c 

0. 

20 

0.23 

0.33 

0.35 

0. 

26 

0 

27 

0. 

35 

0.41 

0 

36 

rionic 

+1. 

03 

+1.53 

+2.64 

-2.93 

+2. 

07 

-5. 

56 

+3 

63 

+3.75 

+3 

61 

relec 

-2. 

73h 

-3.201 

-4.68j 

-4.77k 

-6. 

4/ 

+2. 

66m 

-5 

63k 

-6. 71k 

-7 

40j 

theo 

r 

sum 

-1. 

7 

-1.7 

-2.0 

-1.8 

-4. 

3 

-2. 

9 

-2 

0 

-3.0 

-3 

8 

x* 

exptl 

-1. 

6h 

-l.ln 

2. 0C 

1. 6C 

4. 

3C 

-2. 

4  C 

1. 

8C 

3.0C 

4. 

3C 

NOTE:  a  =  aQ  or  a  r's  represent  r^4  (zincblende)  and  ^3 

(vurtzite)  and  are  in  units  of  10-12  m/v.  rexpt3  are  measure¬ 
ments  with  clamped  crystals.  Their  signs  have  not  been  deter¬ 
mined  unless  specified. 
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Table  A-II.  Comparison  of  the  Geometrical  factors  Between 
Zincblende  and  Wurtzite  and  LiNbO^ 


LiNb03 

Zincblende 

Wurtzite 

(Nb-O)  Short 

(Nb-O)  Long 

£  a  ~ 
n  n3 

0 

0 

2.849 

-4.014 

r  2 

E  a  . 
n  m 

16/3 

8/3 

2.0 

2.0 

Z  a  ,a  0a 
n  al  n2  n3 

16/3  v/T 

0 

0 

0 

2  2 
£  a  ,a  ~  =  £a  0a  « 
n  nl  n3  n2  n3 

0 

-8/9 

1.103 

-1.109 

y  3 
i  a  - 
n  n3 

0 

16/9 

0.642 

-1.797 

y  3 
^  a  0 
n  n2 

0 

0 

0.396 

-0.195 

When 


than  the  electronic  term.  This  is  the  case  in  LiNbO^  and  LiTaO^. 


£  a 

n  m 


0,  as  in  zincblende  and  wurtzite,  we  have  to  settle  for  the 


smaller  term  £  fa  .a  .a  ,  .  This  is  the  main  reason  why  LiTaCL  has  r„~  = 
1  „  n  ni  nj  nk  „  3  33 

30.3  x  10  m/V,  while  GaAs  and  ZnS  have  r,.  -  2  x  10  m/V. 

* 

The  effective  charge  e  is  related  to  the  Szigetti  effective  charge 
*  c 
e  by 
s 

e  +  2  ^ 


The  value  of  e  varies  from  0.2  to  0.4  and  seems  to  be  independent  of 
c 

the  number  of  valence  electrons.  The  calculation  of  r  using  Eqs.  17 

J 

and  18  is  shown  in  Table  I  as  r.  .  .  The  pure  electronic  contribution 

ionic 

is  entered  as  r  ,  .  It  is  obtained  directly  from  the  second  harmonic 

elec  2 

generation  coefficient  by  r^  =  -  4d^/e  .  For  most  of  the  crystals  in 
Table  I,  rionic  positive  and  relec  is  negative.  Therefore,  the  pre¬ 
dicted  electrooptic  coefficients,  rtlieo  =  r.  .  +  r  ,  ,  involve  the 

r  sum  ionic  elec 

algebraic  additi-m  of  two  numbers  of  comparable  values.  The  only 
exception  is  CuCl.  Due  to  the  unfilled  shell  in  Cu,  the  sense  of  bond 
polarization  in  CuCl  is  different  from  that  in  other  crystals.^  For 
example,  excluding  valence  electrons.  As  is  +5,  Ga  is  +3,  and  Cl  is  +7, 
but  Cu  is  +11.  So  the  signs  of  r.  .  and  r  ,  of  CuCl  are  different 
from  those  of  other  crystals.  However,  since  the  magnitude  of  rionic 
is  larger  than  that  of  re^ec  in  CuCl,  we  still  obtain  a  negative 
electrooptic  coefficient  for  CuCl.  Predicted  values  in  Table  I  show 
a  good  agreement  with  experiment.  The  worst  case  is  GaP.  It  is 
interesting  to  note  that  the  electronic  contribution  is  about  double  the 
ionic  contribution.  This  is  in  excellent  agreement  with  the  experimental 
observation. ^ 

We  find  that  the  contributions  to  r„k  from  the  homopolar  part  and 
the  heteropolar  part  are  comparable.  We  intend  to  extend  this  model  to 
complex  crystals  with  different  point  group  symmetries.  The  generaliza¬ 
tion  of  the  bond  parameters  used  above  to  multibond  crystals  has  already 


been  considered.  The  one  parameter  that  will  need  consideration  is  e 


the  effective  ionic  charge.  It  was  found  empirically  to  be  equal  to 


C/fico  in  diatomic  crystals.  If  this  relation  and 
P 


e  +  2 


e  = 
c 


3c 


are  valid  in  the  more  complex  crystals,  then  our  model  can  be  applied  to 
these  cases.  Present  calculations  on  LiNbO^,  LiTaO^,  and  ternary  chal- 
copyrlte  crystals  are  reported  in  this  document  (see  Table  1). 


REFERENCES 


1.  B.F.  Levine,  Phys.  Rev.  Lett.  22,  787  (1969);  25,  440  (1970). 

2.  J.C.  Phillips,  Phys.  Rev.  Lett.  20^,  550  (1968). 

3.  J. A.  Van  Vechten,  Phys.  Rev.  182,  891  (1969);  187,  1007  (1969). 

4.  C.  Kittel,  Introduction  to  Solid  State  Physics,  5th  ed.  (John 

Wiley,  N.Y. ,  1976). 

5.  B.F.  Levine,  Phys.  Rev.  117,  2600  (1973). 

6.  H.B.  Callen,  Phys.  Rev.  76,  1394  (1949). 

7.  I.P.  Kaminow  and  E.H.  Turner,  Phys.  Rev.  B5,  1564  (1972). 

8.  W.L.  Faust  and  C.H.  Henry,  Phys.  Rev.  Lett.  17^,  1265  (1966). 

9.  P.  Lawaetz,  Phys.  Rev.  Lett.  _26,  697  (1971). 


